DGMRES method augmented with eigenvectors for computing the Drazin-inverse solution of singular linear systems

DGMRES: A GMRES-type algorithm for Drazin-inverse solution of singular nonsymmetric linear systems

In a recent work by the author [Linear Algebra Appl. 298 (1999) 99] Krylov subspace methods were derived for the Drazin-inverse solution of consistent or inconsistent linear systems of the form Ax = b, where A ∈ CN×N is a singular and in general non-hermitian matrix that has an arbitrary index. One of these methods, modeled after the generalized conjugate residual method (GCR) and denoted DGCR,...

Orthogonal polynomials and semi-iterative methods for the Drazin-inverse solution of singular linear systems

In this work we present a novel class of semi-iterative methods for theDrazin-inverse solution of singular linear systems,whether consistent or inconsistent. The matrices of these systems are allowed to have arbitrary index and arbitrary spectra in the complex plane. The methods we develop are based on orthogonal polynomials and can all be implemented by 4-term recursion relations independently...

GGMRES: A GMRES--type algorithm for solving singular linear equations with index one

In this paper, an algorithm based on the Drazin generalized conjugate residual (DGMRES) algorithm is proposed for computing the group-inverse solution of singular linear equations with index one. Numerical experiments show that the resulting group-inverse solution is reasonably accurate and its computation time is significantly less than that of group-inverse solution obtained by the DGMRES alg...

A Bi-CG type iterative method for Drazin-inverse solution of singular inconsistent nonsymmetric linear systems of arbitrary index

Convergence of product integration method applied for numerical solution of linear weakly singular Volterra systems

We develop and apply the product integration method to a large class of linear weakly singular Volterra systems. We show that under certain sufficient conditions this method converges. Numerical implementation of the method is illustrated by a benchmark problem originated from heat conduction.